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Abstract 

00 ■ 

Consider in L 2 (M) the Schrodinger operator family H(g) := —d 2 x + V g (x) depending 
on the real parameter <?, where V g (x) is a complex- valued but PT symmetric periodic 
potential. An explicit condition on V is obtained which ensures that the spectrum 
of H(g) is purely real and band shaped; furthermore, a further condition is obtained 
j which ensures that the spectrum contains at least a pair of complex analytic arcs. 

On 

1 Introduction and statement of the results 

There is currently an intense and ever increasing activity on an aspect of quantum theory 
known as PT-symmetric quantum mechanics (see e.g. jBelj . |BBM| . jBe2| . [BBJj . [CJT] . 
[CJN] . [LZ], [Sp], [ZCBR]). Mathematically speaking, in the simplest, one -dimensional 



C3 



in 

o 

oo 



case one deals with the stationary Schrodinger equation 

> 

l>. H*l) :={--—; + V)il) = Ei) . (1.1 



dx 2 

where the potential V(x) can be complex-valued but is invariant under the combined action 
of the linear parity operation P, Pip{x) = ip(—x), and of the (anti)linear "time-reversal" 
symmetry, i.e. the complex-conjugation operation Ttp(x) = namely, V(—x) = 

V{x). The basic mathematical problem is to determine under what conditions, if any, on 
the complex PT symmetric potential V the spectrum of the corresponding Schrodinger 
operator is purely real. 

Here we deal with this problem in the context of periodic potentials on R, already 
considered in [Ahj . [BDMj . fCe) . [CR], [Jo], |Sh] . Without loss, the period is assumed to 
be 2ir. If V is periodic and real valued it is well known (see e.g. [BS| ) that, under mild 
regularity assumptions, the spectrum is absolutely continuous on R and band shaped. 

It is then natural to ask whether or not there exist classes of PT symmetric, complex 
periodic potentials generating Schrodinger operators with real band spectrum. This ques- 
tion has been examined in [Ah], [BDMj . fCe) . |CR| . [Jo], by a combination of numerical 
and WKB techniques, in several particular examples. It was later proved in [Sh] that the 
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above arguments cannot exclude the occurrence of complex spectra, and actually a con- 
dition has been isolated under which H admits complex spectrum consisting of a disjoint 
union of analytic arcs [Sh| . 

Our first result is the explicit determination a class of PT symmetric, complex periodic 
potentials admitting real band spectrum. Denote T(u) the non-negative quadratic form 
in L 2 (R) with domain i? 1 (R) defined by the kinetic energy: 

T(u) := J \u'\ 2 dx, ueH 1 (R) (1.2) 

Let q be a real-valued, tempered distribution. Assume: 

(1) q is a 2-7r-periodic, P symmetric distribution belonging to HT~ (M); 

(2) W(x) :R—>-C belongs to L°°(M) and is PT-symmetric, W(-x) = -W(x), . 

(3) q generates a real quadratic form Q(u) in L 2 (R) with domain H l (M); 

(4) Q(u) is relatively bounded with respect to T(rt) with relative bound b < 1, i.e. there 
are b < 1 and a > such that 

Q(u) < bT(u) +a\\u\\ 2 (1.3) 

Under these assumptions the real quadratic form 

H (u) := T(u) + Q(u) u€H\R) (1.4) 

is closed and bounded below in L 2 (R). We denote H(0) the corresponding self-adjoint 
operator. This is the self-adjoint realization of the formal differential expression (note the 
abuse of notation) 

H(0) = -^+q(x) 
ax z 

Under these circumstances it is known (see e.g. [AGHKH]) that the spectrum of H(0) is 
continuous and band shaped. For n = 1,2,... we denote 

the bands of H(0), and A n :=]P2m (hn+i[> ]®2n+i, ®2(n+2) [ the gaps between the bands. 
Here: 

< «0 < Po < Pi < Oil < «2 < 02 < 03 < «3 < OiA < • • • • 

The maximal multiplication operator by is continuous in L 2 , and therefore so is the 
quadratic form (u,Wu). It follows that the quadratic form family 

H g (u) :=T{u) + Q{u)+g(u,Wu), u G ff 1 (R) (1.5) 
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is closed and sectorial in L 2 (R) for any g G C. 

We denote H(g) the uniquely associated m— sectorial operator in L 2 (R). This is the 
realization of the formal differential operator family 

d 2 

H(g) = + q(x) + gW{x) 

By definition, H (g) is a holomorphic family of operators of type B in the sense of Kato 
for g £ C; by (1) it is also PT symmetric for g G R. Our first result deals with its spectral 
properties. 

Theorem 1.1 Xei a// gaps of H(0) be open, namely: a n < f3 n < a n+ i Vn G N, and /e£ 
i/iere exist d > suc/j i/iaf 

- inf A n := d > (1.6) 

Then, if 

l9l< WTdmu :=v (L7) 

i/iere exist 

< ao(g) < Po(g) < Px{g) < a^g) < a 2 {g) < (3 2 {g) . . . 

sttc/t i/tai 

a(i?( ff )) = ( (J B 2n (g)) |J ( U ^WG/)) (1.8) 

VnGN / VneN / 

where, as above 

E>2n{g) ■= [a2n(s),/?2n(fi0], #2n+l(s0 — [#2n+l (fiO , « 2 n+l (ff)] • 

Remark 

The theorem states that for \g\ small enough the spectrum of the non self-adjoint operator 
H(g) remains real and band-shaped. The proof is critically dependent on the validity of 
the lower bound (|1.6p . Therefore it cannot apply to smooth potentials q(x), in which case 
the gaps vanish as n — > oo. Actually we have the following 
Example 

A locally i? _1 (K) distribution q(x) fulfilling the above conditions is 

q(x) = 5(x — 2im) 

the periodic 5 function. Here we have: 

Q( u ) = \u(2irn)\ 2 = j q(x)\u(x)\ 2 dx, u G H 1 (R) 
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This example is known as the Kronig-Penney model in the one-electron theory of solids. 
Let us verify that condition (II. 3D is satisfied. As is known, this follows from the inequality 
(see e.g. [Ka], §VI.4.10): 

K2vm)| 2 < e / \u'(y)\ 2 dy + 5 / \u(y)\ 2 

J 2irn J 2nn 

where e can be chosen arbitrarity small for 5 large enough. In fact, if u G ff 1 (]R) this 
inequality yields: 

/ q(x)\u(x)\ 2 dx = J2H2™)\ 2 <e [ \u'{y)\ 2 dy + 8 [ \u(y)\ 2 
Jr ngz Jr Jr 

= eT(u) +5\\u\\ 2 , ueH 1 ^) 

which in turn entails the closedness of T{u) + Q(u) defined on H 1 by the standard Kato 
criterion. The closedness and sectoriality of TL g {u) defined on iif 1 (K) is an immediate 
consequence of the continuity of W as a maximal multiplication operator in 1? . For 
the verification of (|1.6|) . see e.g.[AGHKHj. Hence any bounded PT-symmetric periodic 
perturbation of the Kronig-Penney potential has real spectrum for g £ R, |g| < ~g, where 
g is defined by (|1 .T[) . . 

As a second result, we show that an elementary argument of perturbation theory allows 
us to sharpen the result of |Shj about the existence of complex spectra for PT-symmetric 
periodic potentials. 

Let indeed W(x) £ L°°(K;C) be a 27r-periodic function. Then the continuity of W as a 
multiplication operator in L 2 (R) entails that the Schrodinger operator 

d 2 u 

K(g)u:=-—^+gWu, u £ D(K) := H 2 (R), g£R (1.9) 
dx z 

is closed and has non-empty resolvent set. Consider the Fourier expansion of W: 
W{x) = Y w n e inx , w n = — [ W W{x)e' inx dx 

which converges pointwise almost everywhere in [— 7r, it]. Then we have 



Theorem 1.2 Let W(x) be PT '-symmetric, namely W(— x) = W{x). Then 

1. w n = w n , Vn £ Z; 

2. Furthermore, let there exist k £ N, k odd, such that WkW-k < 0. 

Then there is 5 > such that for \g\ < 5 the spectrum of K(g) contains at least a pair of 
complex conjugate analytic arcs. 
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Remarks 

1. This theorem sharpens the results of [Sh] in the sense that its assumptions are 
explicit because they involve only the given potential W(x), while those of Theorem 
3 and Corollary 4 of |Sh] involve some conditions on the Floquet discriminant of the 
equation K(g)ip = Eip. This requires some a priori information on the solutions of 
the equation itself. 

2. Explicit examples of potentials fulfilling the above conditions are: 

W(x) = i sin 2fc+1 nx, k = 0, 1 ... ; n odd 

For g = 1 these potentials have been considered in [ EDM ] ■ where is is claimed that 
the spectrum is purely real. A more careful examination by [Sh| shows that the 
appearance of complex spectra cannot be excluded. 

2 Proof of the statements 

Let us first state an elementary remark under the form of a lemma. Incidentally, this also 
proves Assertion 1 of Theorem 1.2. 

Lemma 2.1 Let f(x) G L°°(R;C) be 2n periodic, f(x + 2vr) = f(x), x G M, and PT 
symmetric, f(—x) = f(x). Consider its Fourier coefficients 

fn = ^- 1^ f(x)e- inx dx 

Then f n = f n Vn G Z. 
Proof 

The assertion is an immediate consequence of the Carleson-Fefferman theorem, which 
states the pointwise convergence of the Fourier expansion 

J(-X) = YJn^ X = E fne mX = f(x) 

almost everywhere in [0, 2ir]. This proves the Lemma. 

To prove Theorem 1.1, let us first recall that by the Floquet-Bloch theory (see e.g. 
[BSJ, [Ea] ) . A G a{H{g)) if and only if the equation H(g)tp = Aip has a non-constant 
bounded solution. In turn, all bounded solutions have the (Bloch) form 

^ p {x;\g) = e i ^<l> p {x-\g) (2.1) 
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where p G] — 1/2,1/2] := B (the Brillouin zone) and <j) p is 27r-periodic. It is indeed 
immediately checked that ip p (x;X,g) solves H(g)ip = Xip if and only if cj) p (x;X,g) is a 
solution of 

Hp(g)4> p (x; A, g) = X<p p (x; X, g). 
Here H p (g) is the operator in L 2 (0, 2-7r) given by 



d 
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H p (g)u= ^-i—+pj u + qu + gWu, u G D(H p (g)) (2.2) 

with periodic boundary conditions; its realization will be recalled below. More precisely, 
denote S 1 the one-dimensional torus, i.e. the interval [—it, it] with the endpoints identified. 
By Assumptions (1) and (2) above the restriction of q to S , still denoted q by a standard 
abuse of notation, belongs to H~ 1 (S 1 ) and generates a real quadratic form Q p (u) in L 2 (5 1 ) 
with domain H 1 ^ 1 ). By Assumption (3) Q p (u) is relatively bounded, with relative bound 
zero, with respect to 

T p (u) := r [-iu ! + pu][iu' + pu] dx, D(T p (u)) = H 1 (S 1 ) (2.3) 

J — 7T 

so that the real semibounded form TL p (u) := T p {u) + Q p {u) defined on H 1 (S 1 ) is closed. 
The corresponding self-adjoint operator in L 2 (5' 1 ) is the self-adjoint realization of the 
formal differential expression (note again the abuse of notation) 

d 2 d 9 

H p (0) = -^-2ip—+p 2 + q 

As above, the form Tt p (g)(u) := T~t p (u) + (u, Wu) defined on i? 1 (5 1 ) is closed and sectorial 
in L 2 (5 : ). Let H p (g) the associated m-sectorial operator in L 2 (5' 1 ). On u G D(H p {g)) the 
action of the operator H p (g) is specified by (|2.2|) : moreover, H p (g) has compact resolvent. 
Let 

°{ h p{9)) ■= {^n(g;p) ■ n = 0, 1, . . .} 

denote the spectrum of H p (g), with p 6] - 1/2, 1/2], \g\ < g. By the above remarks we 
have 

a(H(g))= |J a(H p (g)) = [jX n (g;p) 

pe]-l/2,l/2] P,« 

To prove the reality of a(H(g)), \g\ < g, it is therefore enough to prove the reality of all 
eigenvalues X n {g;p) : n = 0, 1, . . . ;p g] — 1/2, 1/2]; n = 0, 1, — 

To this end, let us further recall the construction of the bands for g = 0: it can 
be proved that under the present conditions all eigenvalues X n (0;p) are simple V p G 
[—1/2, 1/2]; the functions X n (0;p) are continuous and even in [—1/2, 1/2] with respect to 



6 



p, so that one can restrict to p £ [0, 1/2]; the functions A 2 fc(0;p) are strictly increasing on 
[0, 1/2] while the functions A 2 £;+i(0;p) are strictly decreasing, k = 0, 1 Set: 

a fc = A fe (0,0), & = A fe (0, 1/2) 

Then: 

a < A) < #L < «i < «2 < < (k ■ ■ ■ 

The intervals [a 2n , fan] and [/3 2 „+i, ct2n+i] coincide with the range of A 2n (0,p), A 2 „+i(0,p), 
respectively, and represent the bands of o~(H(0)); the intervals 

A„ :=]/32n,/?2n+l[, ]«2n+l j «2(n+2) [ 

the gaps between the bands. 

The monotonicity of the functions A n (0,p) and Assumption (jl.6p entail 

inf min |A„(0,p) - A„+i(0,p)| > 2d (2.4) 

n pe[0,l/2] 

Let us now state the following preliminary result: 
Proposition 2.2 

1. Let g £ T>, where T> is the disk {g : \g\ < g}. For any n, there is a function 
X n (g,p) - T> x[0, 1/2] — * C, holomorphic in g and continuous in p, such that X n (g,p) 
is a simple eigenvalue ofTC p (g) for all (g,p) 6 T> x [0, 1/2]. 

2. 

_sup \K(g,p) - K(o,p)\ < -x 

geX>,p6[0,l/2] 

3. If g € R fl T> all eigenvalues X n (g,p) are real; 

4. IfgGRHV then a{H p {g)) = {X n (g , p)}%> =0 . 

Assuming the validity of this Proposition the proof of Theorem 1.1 is immediate. 
Proof of Theorem 1.1 

Since the functions X n (g;p) are real and continuous for g G R n V, p G [0, 1/2], we can 
define: 

«2n(ff) := min An (5; p); fon(g) ■= max A 2n (#;p) (2.5) 
pe[o,i/2] pe [0,1/2] 

a 2n +i{g) ■= max A 2n +i(ff;p); /Wi(sO := min A 2n+1 (g;p) ( 2 -6) 

p6[0,l/2] pe[0,l/2] 
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Then: 

oo 

a(H(g)) = |J B n {g) 
n=0 

where the bands B n (g) are defined, in analogy with the g = case, by: 

B 2n (g) := [a 2n (g),f32n(g)] 
B2n+i{g) ■= [P2n+i(g),a2n + 1(g)) 

By Assertion 2 of Proposition 2.2 we have, V n = 0, 1, . . ., V <7 £ D: 
a2n(g) - ~ = A 2n (0, 0) - ^ < X 2n (0,p) - ~ 

< X2n(g,p) < A 2 n(0, 1/2) + I = 2n + ~ 



whence 



This yields: 



«2n(s) - g - ^n(g,P) < 02n + Vn, Vp G P. 



«2n — -, P2n + ~ 



d d 

P2n+l - ~, «2n+l + ~ 



B 2 n(g) c 

By an analogous argument: 

B2n+i{g) C 

Therefore the bands are pairwise disjoint, because the gaps 

A n (g) ■=\hn{g),ot2n(g)[ 1 }a 2n +i{g),02n+i(g)[ 

are all open and their width is no smaller than d. In fact, by (jl.6p we have: 

\a n - a n+1 \ > 2d, \(3 n - (3 n+1 \ < 2d. 

This concludes the proof of the Theorem 

We now prove separately the assertions of Proposition 2.2 
Proof of Proposition 2.2, Assertions 1 and 2 

Since the maximal multiplication operator by W is continuous in L 2 (S l ) with norm || W||oo, 
the operator family TC p (g) is a type-A holomorphic family with respect to g £ C, uniformly 
with respect to p £ [0, 1/2]. Hence we can direct apply regular perturbation theory (see 
e.g.[Ka]): the perturbation expansion near any eigenvalue X n (0;p) of 7i(0,p) exists and is 
convergent in T> for g £ T> to a simple eigenvalue X n (g;p) of H(g;p): 



Xn(g;p) = Xn{0,p) + J2\ s n (0;p)g s , g^v 



(2.7) 



s=l 



The convergence radius r n (p) is no smaller than g. Hence g represents a lower bound for 

r n (p) independent of n and p. Moreover A*(0;p) is continuous for all p £ [0,1/2], and 

hence the same is true for the sum X n (g;p). This proves Assertion 1. 

To prove Assertion 2, recall that the coefficients A* (0;p) fulfill the majorization (see [Ka], 

§11.3) 

|A*(0;p)| < (——. fJVnn) ^ (\\W\Ud) s (2.8) 

Therefore, by ([2TTD : 

IX ( \ wn M < \9\\\W\\oq 

\K{g;p) - K{o,p)\ < - — , , ... ,, — -ttt = — rnrfiTii - < o 

1 - { 2 \9\\\W Woo/ d) d-lglWWWoo 2 

whence the stated majorization on account of (|1.7|) . 
Proof of Proposition 2.2, Assertion 3 

As is known, and anyway very easy to verify, the PT symmetry entails that the eigenvalues 
of a PT-symmetric operator are either real or complex conjugate. By standard regular 
perturbation theory (see e.g. |Ka| . §VII.2) any eigenvalue A n (0;p) of H p (0) is stable with 
respect to H p (g); since A n (0,p) is simple, for g suitably small there is one and only one 
eigenvalue X n (g,p) of H p (g) near A„(0,p), and X n (g,p) — ► A n (0,p) as g — > 0. This excludes 
the existence of the complex conjugate eigenvalue X n (g,p) distinct from X n (g,p). Thus for 
\g\ suitably small, X n (g,p) is real. This entails the reality of series expansion (12. 7\i 
for g small and hence Vg G P. This in turn implies the reality of X n {g,p)\/ g £ T>. 
Proof of Proposition 2.2, Assertion 4 

We repeat here the argument introduced in [CGS],[CG] to prove the analogous result in 
different contexts. We describe all details to make the paper self contained. We have 
seen that for any rGN the Rayleigh-Schrbdinger perturbation expansion associated with 
the eigenvalue X r (g;p) of H p (g) which converges to X r (0;p) as g — > 0, has radius of 
convergence no smaller than g. Hence, Vg € R such that \g\ < g, H p (g) admits a sequence 
of real eigenvalues X r (g;p),r £ N. We want to prove that for \g\ < ~g, g £ R, H p (g) has 
no other eigenvalues. Thus all its eigenvalues are real. To this end, for any r £ N let Q r 
denote the circle centered at A r (0;p) with radius d. Then if g £ R, \g\ < g, and X(g) is an 
eigenvalue of H p (g): 

Kg) G IJ Sr- 

rGN 

In fact, denoting 

R (z) := (flp(0) - z)- 1 
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for any z £ [J Q r we have 

rGN 

\\gWR (z)\\ < \g\\\W\U\R (z)\\ < gWWlUdistiz^iHo))}- 1 < g|l ^ l|o ° < 1. (2.9) 

The last inequality in (|2.9p follows directly from the definition (|1.7p of 5. Thus, z G 
P( H p{d)) and 

flfo, z) := (H p (g) - zy 1 = R (z)[l + gWRo(z)}- 1 . 

Now let go E K be fixed with \g\ < g. Without loss of generality we assume that go > 0. Let 
X(go) be a given eigenvalue of H p (go). Then X(go) must be contained in the interior (and 
not on the boundary) of Q m for some tiq £ N. Moreover if ttiq is the multiplicity of X(go), 

for g close to 50 there are mo eigenvalues (counting multiplicities) X^(g),s = 1 ,mo, 

of H p (g) which converge to A(#o) as 5 — > go and each function A( s )(g) represents a branch 
of one or several holomorphic functions which have at most algebraic singularities at 
9 = 9o ( see [KaJ, Thm. VII. 1.8). Let us now consider any one of such branches A" (g) 
for < g < go, suppressing the index s from now on. First of all we notice that, by 
continuity, X(g) cannot lie outside Q no for g close to go- Moreover, if we denote T t the 
boundary of the circle centered at X no (0;p) with radius t, < t < d, we have, for z £ Tt 
and < g < g , 

\\gWRo(z)\\ < gllWllooidist&v&pQ)))]- 1 < oW^/i. (2.10) 

Then t > W||oo implies z £ a(H p (g)), i.e. if z e a(H p (g)) n T t then t < g||W||oo < 

50 1 1 W 1 1 00 < 5||W||oo < d. Hence we observe that as g — > ^q", X(g) is contained in the circle 
centered at X no (0;p) and radius g[|W[|oo- Suppose that the holomorphic function X(g) is 
defined on the interval ]<7i,5o] with g\ > 0. We will show that it can be continued up to 
g = 0, and in fact up to g = —~g. From what has been established so far the function 
X(g) is bounded as g —> gf . Thus, by the well known properties on the stability of the 
eigenvalues of the analytic families of operators, X(g) must converge to an eigenvalue X(gi) 
of Hp(gi) as g — > gf and X(gi) is contained in the circle centered at A no (0;p) and radius 

51 Halloo- Repeating the argument starting now from X(gi), we can continue X(g) to a 
holomorphic function on an interval ]g2,gi], which has at most an algebraic singularity at 
g = <?2- We build in this way a sequence g± > gi > . . . > g n > . . . which can accumulate 
only at g = —g. In particular the function X(g) is piecewise holomorphic on ] — g~,~g\. 
But while passing through g = 0, A(g) coincides with the eigenvalue X r (g;p) generated by 
an unperturbed eigenvalue X r (0;p) of H p (0) (namely A„ (0;p)), which represents a real 
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analytic function defined for g g] — g,g[. Thus, A (go) arises from this function and is 
therefore real. This concludes the proof of Assertion 4. 

Proof of Theorem 1.2 

Consider the operator K p (g) acting in L 2 {S l ), defined on the domain H 2 ^ 1 ). By the 
Floquet-Bloch theory recalled above, we have again 

a(K(g)) = |J a(K p (g)). 

pS]-l/2,l/2] 

It is then enough to prove that there is rj > such that K p {g) has complex eigenvalues for 
p G]l/2 — rj, 1/2]. Since K p {g) is PT-symmetric, eigenvalues may occur only in complex- 
conjugate pairs. The eigenvalues of K p (0) are X n (0,p) = (n + p) 2 : n G Z. The eigenvalue 
Ao(0,p) = p 2 is simple Vp G [0, 1/2]; any other eigenvalues is simple for p ^ 0, p ^ 1/2 and 
has multiplicity 2 for p = or p = 1/2 because n 2 = (— n) 2 and (ra+1/2) 2 = (— n— 1+1/2) 2 , 
n = 0, 1, . . .. For p = 1/2 a set of orthonormal eigenfunctions corresponding to the double 
eigenvalue (n + 1/2) 2 = (— n — 1 + 1/2) 2 , n = 0, 1, ... is given by where: 

u n := — ^=e mx , n G Z 

\/27r 

Remark that, for < p < 1/2, u n and U- n -i are the (normalized) eigenfunctions cor- 
responding to the simple eigenvalues (n + p) 2 , (— n — 1 +p) 2 . Let now n = k G N, and 
consider the 2x2 matrix 



T : = 

A trivial computation yields 

T : 



{u k ,Wu k ) (u k ,Wu- k -i) 
(u-k-i,Wu k ) (u- k -x,Wu- k -i) 



w 2k +i 

W-2k-l 



with purely imaginary eigenvalues ^ = iziy/— W2 k +i'W-2k-i- By standard degenerate 
perturbation theory, for g small enough TCi/2(g) admits the pair of complex conjugate 
eigenvalues 

Xf(g, 1/2) = (k + 1/2) 2 ± i gy /-w2k+iw-2k-i + 0(g 2 ) (2.11) 

Under the above assumptions, for any fixed g£K the operator family p ^ K p (g) is type- A 
holomorphic in the sense of Kato (see |Ka| . §VII.2) \/ p G C because its domain does not 
depend on p and the scalar products (u, K p (g)u) : u G H 1 ^ 1 ) are obviously holomorphic 
functions of p \/p G C. This entails the continuity with respect to p of the eigenvalues 
X n (g;p). Therefore, for \g\ suitably small, there is rj(g) > such that 

lmXi t (g,p)^0, l/2-77<p<l/2 
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It follows (see e.g. [55]. |Eaj ) that the complex arcs £ k := Range(A fc (g,p)) ■ p £ [1/2 — 
77(5), 1/2] belong to the spectrum of K(g). This concludes the proof of the theorem. 
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